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ABSTRACT 

We have made the first direct symmetry tests in the decays of polarized bosons into 
fully-identified bbg states, collected in the SLD experiment at SLAC. We searched for evi- 
dence of parity violation at the bbg vertex by studying the asymmetries in the 6-quark polar- 
and azimuthal-angle distributions, and for evidence of T-odd, CP-even or odd, final-state 
interactions by measuring angular correlations between the three-jet plane and the po- 
larization. We found results consistent with Standard Model expectations and set 95% C.L. 
limits on anomalous contributions. 
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Polarization is an essential tool in investigations of fundamental symmetries in 
elementary-particle interactions. Parity violation was first discovered using jS decays 
of polarized ^°Co and T, CP and CPT violations were searched for using polar- 
ized neutrons [0 and polarized positronium 0. The unique sample of polarized 
bosons produced in annihilations of longitudinally-polarized electrons with unpolarized 
positrons at the SLAG Linear Collider (SLC) can similarly be employed for fundamen- 
tal symmetry tests. Here we use e"'"e~ ^ —>■ three-jet events to test symmetry 
properties of the Standard Model (SM). The bbg final state provides a particularly 
interesting probe for possible beyond-SM processes that couple to massive particles. 

The tree-level differential cross section for e"'"e~ —>■ qqg can be expressed ^ 
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where 6 is the polar angle of the thrust axis w.r.t. the electron beam, and x is the 
azimuthal angle of the event plane w.r.t. the quark-electron plane. The thrust axis is 
defined to be along the most energetic jet and to point into (opposite) the hemisphere 
containing the quark (antiquark) if the quark (antiquark) has the higher energy. The 
sign of cosx is defined in terms of momentum vectors, sign(cosx) = sign((gx^)-(gxe^)). 
The functions h^j^ and h^p contain the dependence on the beam polarization and the 
electroweak couplings [Q; QCD contributions are expressed in terms of cxj = d'^ai/dxdx, 
with i = U,L, T, I, A, P, where x and x are the scaled momenta of the quark and anti- 
quark, respectively. While the first four terms are even under P reversal, the last two 
terms are P-odd, and are sensitive to any parity-violating interactions at the Z^qq or 
qqg vertices. More generally there can be, in addition, three terms that are odd under 
T reversal . However, these terms vanish at tree level in a theory that respects CPT 
invariance. 

We report the first experimental study of angular asymmetries in polarized Z^ 
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decays to fully-identified bbg states. We present new tests of QCD using the P-odd 
asymmetries in the 6 and x distributions as well as the first measurements of 
two T-odd triple-product correlations between the bbg plane and the polarization. 
Integrating Eq. |l] over x, x and Xi 

— ^ ocl + a cos^ e + 2Pz Ap cos 9, (2) 
a cosy 

where Pz = {Pe- — ^e) /(I ~ Pe- ■ ^e); -Pg" IS the signed electron-beam polarization, the 
parity- violation parameter Ap = Af,-crp/{(ru + crL), A^ (Af,) is the electroweak coupling 
of the to the initial (final) state, and Aj = 2vjaj/{vj + a|) in terms of the vector, 
Vj, and axial-vector, aj, couplings of fermion j to the Z^; a = {au — 2d-L)/{au + ^ai), 
and (Ti = j d?ai. Similarly, 

d(7 Stt 

— oc 1 + /3cos2x - ^-^ Pz A'p cosx, (3) 

where A'p = Ah ■ o-a/^o-u + and [3 = crT/{d'u + ctl)- Given the value of A},, 
measurement of Ap and A'p allows one to test the QCD predictions for opjiau + di) 
and oaH^^u + Furthermore, the ratio J^pjAp yields oaIop independently of A^,. 

In terms of the polar angle cu, w.r.t. the electron-beam direction, of the vector n 
normal to the event plane, 

da ^ . 2 16 



oc 1 -|- 7 cos uj ^ — Pz Atcoscj, (4) 
d cos ijj 9 

where 7 = {2a l — cfu — 6o-t)/(3o-(/ + 2di^ + 2aT)\ the asjTiimetry term is one of the 
three T-odd terms mentioned above. The vector n can be defined in several ways; for 
example: 1) using the two highest-energy jets, n = pi x or 2) using the quark and 
anti-quark momenta, n = ph x p^. The asymmetry is CP-even in the first definition, 
and CP-odd in the second; in both cases in the SM At = at tree level. Higher-order 
corrections to e^e~ -^hbg yield \A'j'\ < 10^^ 0. The asymmetry in costu is hence 
potentially sensitive to beyond-SM processes [Q. 

The measurement was performed using approximately 550,000 Z^ hadrons 
decays produced in collisions of longitudinally-polarized electrons with unpolarized 



positrons at the SLC between 1993 and 1998. In order to reduce systematic effects on 
polarization-dependent asymmetries the electron polarization direction was reversed 
randomly pulse-by-pulse; the magnitude of the average polarization was 0.73. The data 
were recorded in the SLC Large Detector (SLD) 0. The trigger and hadronic event 
selection criteria are described elsewhere [jlOl- This analysis used charged tracks mea- 
sured in the central drift chamber and in the CCD-based vertex detectors (VXD) ||11|| . 
About 70% of the data were taken with the new VXD installed in 1996 (VXD3), and 
the rest with the previous detector (VXD2). Only well-reconstructed tracks |]T2[ were 
used for the 6-jet tagging. Particle energies were measured in the liquid-argon and 
warm-iron calorimeters; the calorimetric information was used only to reconstruct the 
thrust axis. 



In each event jets were reconstructed using the "Durham" algorithm ||13[. To select 
planar three-jet events we required exactly three reconstructed jets to be found with a 
jet-resolution parameter value of ?/c=0.005, the sum of the angles between the three jets 
to be greater than 358°, and that each jet contain at least two charged tracks; 74,886 
events satisfied these criteria. The jet energies were calculated by using the measured 
jet directions and solving the three-body kinematics assuming massless jets. The jets 
were then labeled such that Ei > E2 > E^. 

To select bbg events the long lifetime and large invariant mass of 5-hadrons were 
exploited. An algorithm [TB] was applied to the set of well reconstructed tracks in 
each jet in an attempt to reconstruct a decay vertex. Vertices were required to contain 
at least two tracks, and to be separated from the interaction point (IP) by at least 1 
mm. We calculated (described below) that the probability for reconstructing at least 
one such vertex was ~ 91% (77%) in bbg events, ~ 45% (26%) in ccg events, and ~ 2% 
(2%) in light-quark events recorded in VXD3 (VXD2). 

An event was selected as bbg if at least one jet contained a vertex with invariant 
TsU > 1.5 GeV/c^. A total of 14,658 events satisfied this requirement and were 
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subjected to further analysis. We calculated that this selection is 84% (69%) efficient 
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for identifying a sample of bbg events with 84% (87%) purity, and containing 14% 
(11%) ccg and 2% (2%) light-flavor backgrounds. 

In order to isolate the gluon jet, a 6-tag was evaluated for each jet in the selected 
event sample. This tag was defined to be positive if the jet contained a vertex, or if it 
contained at least 3 "significant" tracks, i.e. with normalized impact parameter w.r.t 
the IP d/ad > 3 [|l^]; otherwise the tag was defined to be negative. Jet 1 was tagged 
as the gluon jet (g) if it had a negative 6-tag and both jets 2 and 3 had positive 6-tags. 
Jet 2 was tagged as the gluon-jet if it had a negative 6-tag and jet 3 a positive one. 
Otherwise, jet 3 was tagged as the gluon jet. We calculated that, on average, the purity 
of the tagged gluon jet sample was 91% (88%). 

In order to distinguish the b and b jets we evaluated for each jet j the momentum- 
weighted charge, Qj = J^QilPi ' t\'^, where k=0.5, t is the unit vector along the thrust 
axis and qi and pi are the charge and momentum of the i*^ track in jet j. We calcu- 
lated the charge difference, Qdiff = Qi — Q2 — Qs- If Qdiff was negative (positive) we 
assigned jet 1 as the b (b) jet, unless jet 1 was gluon tagged, in which case we assigned 
jet 2 as the b (b) jet if Qdiff was positive (negative). The probability of correct b/b 
assignment, P(Q[){{), was calculated from the data by using a self-calibration tech- 
nique 0, V{Qi){{) = 00/(00 + ]-°Ll2r>{{l)^ where ab is a fitted parameter. Averaged 
over cos^ at = 0.218±0.018 (VXD3) (0.248±0.030 (VXD2)), corresponding to < V > 
= 0.68 (VXD3) (0.67 (VXD2)). 

A JETSET 7.4 [|l^-based Monte Carlo simulation of hadronic decays combined 
with a simulation of the detector response was used to evaluate the efficiency and purity 
of the bbg event selection and the purity of the jet flavor tags. For those simulated events 
satisfying the three-jet criteria, exactly three jets were reconstructed at the parton level 
by applying the jet algorithm to the parton four-momenta. The three parton-level jets 
were associated with the three detector-level jets by choosing the combination that 
minimized the sum of the angular differences between the corresponding jet axes, and 
the energies and charges of the matching jets were compared. 



5 



Figs. 1(a) and (b) show the observed cos 6* distributions for event samples pro- 
duced by left- and right-handed electron beams, respectively. The distributions may 
be described by 

da 



ocl+acos^^ +2Pz cose [Ap fb {2 pI-1) + 
Ap,, U {2pl - 1) + Ap^uds (1 - A - fc) i2pf^ - 1) ], (5) 



where and fc are the fractions of bbg and ccg events in the sample, respectively, 
Ap^c and Ap^uds are the asymmetry parameters for the respective backgrounds (Fig. 1), 
and pg, pg and pg'^'^ are the probabilities to sign cos6 correctly in bbg, ccg and light- 
quark events, respectively. The data were used to calculate Pg{ = V above); all other 
quantities were calculated from the simulation. A maximum-likelihood fit of Eq. |^ 
yielded Ap = 0.855 ± 0.050 (stat.). 

Figs. 2(a) and (b) show the observed x distributions for events produced with left- 
and right-handed electron beams, respectively. They may be described by 

^ oc 1 + f3cos2x -l^Pz cosx [A^ (2p^ - 1) + 

A'p^, fc {2pl - 1) + (1 - /, - /,) (2pf - 1) ], (6) 



where A'p^ and A'p^^^ are the asymmetry parameters for the respective backgrounds, 
and p^, p'^ and p^^^ are the probabilities to sign cosx correctly in bbg, ccg and light- 
quark events, respectively. Averaged over x, p^ = 0.64 (0.63), which was derived 
using the measured value of V (above) combined with the simulated probability to tag 
the gluon jet correctly. All other parameters were calculated from the simulation. A 
maximum-likelihood fit of Eq. ^ yielded A'p = —0.013 ± 0.033 (stat.). 

Figs. 3(a) and (b) show the left-right-forward-backward asymmetry in | cosu;| = z, 

ApBiz) = ^^^^^ ~ ^Lj-z) + apj-z) - apiz) 
aiiz) + aii-z) + api-z) + apiz) 

for the two definitions of the direction of n: (1) pi x p^, and (2) Ph^pi- No asymmetry 
is apparent. The coso; distributions may be described, assuming no asymmetries in the 
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ccg and light-quark backgrounds, by 



oc 1 cos^ uj H Pz A^fb{2p\ — 1) coso;, (7) 



d cos u; 3 9 

where p\. is the probabihty to sign coscj correctly in the CP-even (+) and CP-odd 
(— ) cases, respectively. In the CP-even case, in which the jets were labeled according 
to their energy, six detector-jet energy orderings were possible for a given parton-jet 
energy ordering. Using the simulation we calculated that, averaged over coscj, p'^ = 
0.76 (0.76). In the CP-odd case, both the gluon jet and the 6-jet must be tagged 
correctly, and p'L ( = above) = 0.64 (0.63). Maximum-likelihood fits (Fig. 3) of 
Eq. yielded = -0.014 ± 0.016 (stat.) and = -0.035 ± 0.024 (stat.). 



We considered sources of systematic error [|1^, [T^ which potentially affect our 
results. The modeling of the detector response was studied by varying the track- 
ing efficiency and resolution within their estimated uncertainties, and by varying the 
vertex-mass and significant-track requirements. In addition, the error on the proba- 
bility for correct jet-energy ordering was estimated from the difference between results 



derived using HERWIG [T^ and JETSET. The error on the 6/6-jet identification prob- 
ability was derived from the statistical uncertainty in the at determination using the 
self-calibration technique. Contributions from the modeling of underlying physics pro- 
cesses were also studied. In bb events we considered the uncertainties on: the branching 
fraction for ^bb, the S-hadron fragmentation function, the rates of production of 
-B^, 5° and mesons, and B baryons, the lifetimes of B mesons and baryons, and 
the average i?-hadron decay charge multiplicity. In cc events we considered the un- 
certainties on: the branching fraction for ^cc, the charmed hadron fragmentation 
function, the rates of production of D^, and Dg mesons, and charmed baryons, 
and the charged multiplicity of charmed hadron decays. We also considered the rate 
of production of ss in the jet fragmentation process, and the production of secondary 
bb and cc from gluon splitting. The uncertainty on the beam polarization was also 
taken into account. 

The systematic errors on A'p and are negligibly small as the uncertainty dimin- 



ishes with the measured asymmetry, which in all three cases is small and consistent with 
zero. In the case of the largest error contributions arose from the ah uncertainty 
{SApjAp ~5.7%) and from the uncertainty in the background level {SApjAp ~3.9%); 
other contributions were at the level of 1% or less. 

From the measured values of and A'p, and assuming the SM expectation of A^ 
^ 0.935 for sin^ 6*^=0. 232, we find, respectively 

, = 0.914 ± 0.053 {stat.) ± 0.063 {syst.), 

cru + (Jl 

^ ^ = -0.014 ± 0.035 (stat.) ± 0.002 (syst). 

These values are consistent with the Ola^) QCD expectation (for massless quarks) 
of d-p/{au + = 0.93 and aAl^o-u + = —0.06, respectively, calculated using 
JETSET 7.4. These yield (JaIop = —0.015 ± 0.038 independent of the assumed A^,, 
and consistent with the expected value —0.065. We also used the Ap and A'p values 
to set a 95% CL. limit on an anomalous axial-vector coupling of the gluon to the 
6-quark, parametrized by a factor (1 + 675)7,^ in the bbg coupling, of e < 0.34. The 
measured values of Aj^ and A^ correspond to 95% CL. limits of —0.045 < Aj. < 0.016 
and -0.082 < A^ < 0.012 (Fig. 3). 

In conclusion, we have made the first symmetry tests in polarized Z° decays to 
bbg. From the forward-backward polar-angle asymmetry of the signed-thrust axis and 
the azimuthal-angle asymmetry we found the parity- violation parameters Ap and Ap, 
respectively, to be consistent with 0(«^) QCD expectations. We set a corresponding 
95% CL. limit on parity violation at the bbg vertex. Using the event-plane normal 
polar-angle distributions we set 95% CL. limits on the T-odd, CP-even and odd asym- 
metry parameters Aj, and A^, respectively. 

We thank Arnd Brandenburg, Lance Dixon and Per Osland for many helpful discus- 
sions. We thank the personnel of the SLAC accelerator department and the technical 
staffs of our collaborating institutions for their outstanding efforts on our behalf. 
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Figure 1. Polar-angle distribution of the signed-thrust axis direction (see text) for 
(a) left- and (b) right-handed electron beam. The histograms represent the simulated 
backgrounds. 
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Figure 2. As Fig. 1, for the azimuthal-angle distribution (see text). 
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Figure 3. Left-right forward-backward asymmetry in cosu (see text) for (a) CP-even 
case, and (b) CP-odd case. The sohd curve is the best fit to the data sample, and the 
dashed curves correspond to the 95% C.L. hmits. 
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